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ABSTRACT 

We evaluate the ability of Spider, a balloon-borne polarimeter, to detect a divergence-free polariza- 
tion pattern ( U B- modes" ) in the Cosmic Microwave Background (CMB). In the inflationary scenario, 
the amplitude of this signal is proportional to that of the primordial scalar perturbations through 
the tensor-to-scalar ratio r. We show that the expected level of systematic error in the Spider in- 
strument is significantly below the amplitude of an interesting cosmological signal with r — 0.03. 
We present a scanning strategy that enables us to minimize uncertainty in the reconstruction of the 
Stokes parameters used to characterize the CMB, while accessing a relatively wide range of angular 
scales. Evaluating the amplitude of the polarized Galactic emission in the Spider field, we conclude 
that the polarized emission from interstellar dust is as bright or brighter than the cosmological signal 
at all Spider frequencies (90 GHz, 150 GHz, and 280 GHz), a situation similar to that found in the 
"Southern Hole." We show that two ~ 20-day flights of the Spider instrument can constrain the 
amplitude of the i?-mode signal to r < 0.03 (99% CL) even when foreground contamination is taken 
into account. In the absence of foregrounds, the same limit can be reached after one 20-day flight. 
Subject headings: balloons - cosmic background radiation - instrumentation: polarimeters 



1. INTRODUCTION 

Over the past decade, our understanding of cosmology 
has been revolutionized by ever more precise measure- 
ments of the temperature and of the polarization of the 
Cosmic Microwave Background (CMB). The first-year 
WMAP data provided the first full-sky characterization 
of the C MB temperature an isotropics, first detected by 
COBE (Smoot et al. 1992) and mapped with high fi- 
delity by BOOMERANG Qde Bernardis et al.|2000[ ), over a 
range of angular scales that allowe d tor the derivation o f 
stringent cosmological constraints ( |Spergel et al. 2003). 
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These constraints were further strengthened by WMAP's 
full-sky measurement of a c url-free polarizat ion pattern 
(£-mode) in the CMB sky ( |Page et al ||2007D , a pattern 
first detected from the ground by DAS1 (Kovac et al. 
2002|). Many other CMB experiments, both ground- 



based and balloon-borne, have provided complementary 
measurements, which made it possible to study physics 
poorly constrained by WM AP alone. T his is the case of 
QUaD ( |Brown et al.|2009[ ) and ACBAR ( |Reichardt et al 



20091, whose data combined with WMAP's enabled the 



convincing detection of the effect of prim ordial helium on 
the C MB temperature power spectrum (Komatsu et al. 



20111, thereby openi ng a new windo w on big b ang nu- 
cleosynthes is. ACT (IDas et al.||2011|) and SPT (IKeisler 



et al have provided temperature data down to ar 



emmute scales, probing the CMB power s pectrum up to 
the s eventh acoustic peak, while BiCEP (Chiang et al. 
20101 first detected the £ ~ 100 acoustic peak in the 
i/'-mode power spectrum, a detection later confirme d by 
QUIET dQuiET Collaboration: Bischoff et al.||2010|). 



Remarkably, the six parameters of the simple flat, 
power-law ACDM model are sufficient to fit not only all 
the available CMB temperature and polarization data, 
but also a wealth of other cosmological probes, such as 
distance mea surements from Bar yon Acoustic Oscilla- 
tions (BAOs, [Percival et al.||2010|, and Hubble co nstant 
measurements from Cepheids (Riess et al. 2009). Al- 
though it is possible to fit more paramet ers to the data, 
there is no indication that any are needed ( Komatsu et al. 



20111. The inflation paradigm, which postulates the ex- 



istence of a period of accelerated expansion in the early 
Universe, is the corner stone of this very successful model 
( Baumann et al.|2009 ) . Not only does it explain the flat- 
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ness and the homogeneity of the Universe, it also provides 
a natural mechanism to generate the primordial Gaussian 
density fluctuations that left their imprint in the form 
of anisotropies in the CMB temperature and in the dis- 
tribution of the large-scale structure we observe today. 
However, although inflation is consistent with current 
observational data, there is no direct evidence to date in 
favor of this paradigm over its viable alt ernatives. The 
ekpyrotic scenario (for a review, see, e.g., |Lehners 2008) 



is arguably the best motivated alternative to inflation. 
In addition to also being consistent with current data, 
it has the added benefit of naturally avoiding some of 
the theoretical complications that one encounters when 
building an inflationary model. 

The predictions of the ekpyrotic scenario and of the in- 
flationary paradigm differ in one major way. The latter 
generically predicts the existence of gravity wave per- 
turbations in the early Universe, which would result in 
an observable divergence-free polarization pattern in the 
CMB (B-modes), whereas the former excludes th e pro 



duction of such perturbations at a detectable level ( Boyle 
et al. 2004). B-modes therefore constitute an unambigu- 
ous signature of inflation. Moreover, the energy scale of 
inflation sets the amplitude of these tensor (gravity wave) 
perturbations relative to that of the scalar (density) per- 
turbations in the form of the tensor-to-scalar ratio r. In 
addition to validating the inflationary paradigm, detect- 
ing B-modes would therefore determine the energy at 
which inflation occurs, which would remove the largest 
source of uncertainty in inflationary model build ing 
B-modes remain as yet undetected. Bicep (Chiang 



et al.| 2010) derived the strongest constraint on r di- 
rectly from B-mode nondetection to date: r < 0.72 
(95% CL). Using CMB temperature data from SPT, 
along wit h CMB temperature a nd polarization data from 
WMAP flKomatsu et al.||2011|), distance measurements 



from BAOs (Percival et al. 



on the Hubble constant (Kiess et al. 2009) 



mop, 



and recent constraints 



(2011 ) found an upper limit of r < 



Keisler et al. 
T7T95% CL). This 
constraint is about an order of magnitude above that nec- 
essary t o rule out the simplest viable class of inflationary 
models ( |Baumann et al.|[2009[ ). 

The Planck satellite has already revolutionized our un- 
derstanding of the sky at microwave to submillim eter 
wavelengths [ Planck Collaboration: Ade et al.||2011 ). It 
has so far proven to be a remarkable machine to study 
the properties of the interstellar medium, and in particu- 
lar of the emission from Galactic dust, which is expected 
to be the brightest polarized sky signal in all channels of 
Planck's High Frequency Instrument (HFI). Even though 
Planck's sensitivity will be sufficient to map the diffuse 
Galactic dust emission and constrain the amplitude of its 
polarized component, it will not be high enough to pro- 
duce high signal-to-noise maps of this polarized emission 
at scales smaller than ~ 20 degrees. This in turn limits 
Planck's ability to probe the cosmological B-mode signal 
beyond the large-scale reionization bump located below 
I ~ 10. Optimistic forecasts of Planck's ability to de- 



tect or constrain .B-mode s have been published (e.g., Ef- 
stathiou fc Gratton|2009 1 . However, two major assump- 



tions are made m these studies: (1) that the Planck HFI 
is ideal at large angular scales (1// noise is for example 
ignored), and (2) that residuals after foreground removal 
are small enough that the induced bias on the determina- 



tion of r is negligible. Although it is difficult to evaluate 
whether this latter assumption is reasonable, the former 
is most li kely problematic (for an early review of th e HFI 
noise, see \Planck HFI Core Team: Ade et al.||2011[ ). 



Spider ( jFilippim et al.||20iup is a balloon-borne po- 
larimeter designed to probe the polarization of the 
microwave sky with unprecedented sensitivity and fi- 
delity. The Antarctic Long Duration stratospheric Bal- 
loon (LDB) platform provides a combination of long 
flight times, near-space optical backgrounds, and sensi- 
tivity to angular scales that are inaccessible from even 
the most favorable ground-based sites. By taking ad- 
vantage of these observing conditions during two flights 
to be launched, respectively, in December 2012 and De- 
cember 2014, Spider will produce high signal-to-noise 
polarized maps of 10% of the sky that are exceptionally 
clean of polarized Galactic emissions, with the goal of 
detecting and characterizing the i ~ 80 acoustic peak in 
the B-mode power spectrum for values of the tensor-to- 
scalar ratio r down to 0.03. In this paper, we review the 
current status of the Spider experiment, and present our 
observing plans to achieve this science goal. 

This paper is organized as follows. In Section [2] we 
briefly describe the Spider instrument, and evaluate the 
level of B-mode contamination exp ecte d from known 
sources of systematic error. Section 2.3 presents a de- 
tailed study of the impact of polarized sidelobes on our 
observations. The expected level of polarized Galactic 
foregrounds in the Spider field is summarized in Sec- 
tion [3j where we also present a model for the polarized 
emission from interstellar dust. In Section [4j we provide 
an overview of our observing strategy, along with a com- 
putation of Spider's filter transfer function Fi. Section[5] 
details the reasoning that led to the selection of 90 GHz, 
150 GHz, and 280 GHz as Spider's observing frequen- 
cies, and gives our expected constraints on r taking into 
account polarized foreground contamination. Our con- 
clusions are summarized in Section [6J 



2. THE Spider INSTRUMENT 
2.1. Description 

The Spider payload (Figure I]) consists of six mono- 
chromatic refracting telescopes derived from the design of 
the highly successful Bicep instrument. This compact, 
axisymmetric optical system greatly limits the opportu- 
nity for polarized systematics. The relatively small aper- 
ture (~ 30' beam FWHM at 150 GHz) allows relatively 
easy baffling of optical sidelobes while retaining sensitiv- 
ity to the angular scales relevant to inflationary science. 
The far-field regime of the optics is also relatively nearby 
(~ 100 m), greatly simplifying optical characterization in 
the laboratory. 

Each refractor illuminates a focal plane of large- 
form at antenna-coup led transition-edge sensor (TES) ar- 
rays (Kuo et al. 2008 ). Each of a focal plane's optical pix- 
els consists of two interpenetrating arrays of slot anten- 
nas sensitive to perpendicular polarization components. 
The power received by each antenna array is coupled 
to a dedicated TES through an inline band-defining mi- 
crostrip filter. The antennas, filters and bolometers for 
hundreds of polarimetric pixels are fabricated photolitho- 
graphically onto a single silicon tile, with no need for ex- 
ternal feed horns or band-defining filters. The arrays are 
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Fig. 1. — Left: Rendering of the Spider payload, showing the cryostat, surrounding frame, and sun shield. Center: Cross-section of 
the Spider telescope CAD model with key components labeled. Top right: Optical pixel, showing phased array of slot antennas feeding 
two TESs at bottom. Bottom right: Close-up of a single TES assembly, showing the TES island and the SiN thermal isolation legs. The 
meandered leg design allows for low thermal conductance in a narrow geometry. 



read out using a time-domain SQUID multiplex er sys- 
tem ride Korte et al.||2003| [Battistelli et al.||2008[ ). Four 
such focal planes are currently in operatio n at the South 
Pole, as part of the Bicep2 instrument riOgburn et al. 



2010| ) and the Keck array ( |Sheehy et aT|2010p 
Spider will deploy telescopes in three treque 



deploy telescopes in three frequency bands 
centered at 90, 150 and 280 GHz. Table [l] lists some of 
the major characteristics of a single Spider receiver in 
each of these observation bands. These frequencies are 
chosen for their sensitivity to the CMB and to Galac- 
tic dust emission. Since all frequency-specific optical el- 
ements are fully contained within each telescope insert, 
Spider's frequency coverage is easily adjustable by swap- 
ping out one or more of these telescopes. This modularity 
gives Spider a great deal of flexibility in choosing an op- 
timal frequency coverage for each flight. Band selection 
is discussed further in Section [H 

The six receivers are housed within a shared ~ 1300 L 
liquid helium cryostat, with sub-Kelvin cooling for each 
focal plane provided by a dedicated 3 He sorption refriger- 
ator. The cryostat is supported below the balloon within 
a lightweight carbon-fiber gondola frame, derived from 
the BLAST ( |Pascale et al.||2008 ) design. A reaction 
wheel allows the payload to scan m azimuth, while the 
elevation of the inner frame is adjustable wi th a simple 



linear actuator similar to the Boomerang (Masi et al. 
2006) design. Pointing information is provided by two 
tracking star cameras, one fixed star camera, rate gy- 
roscopes, differential GPS, and a sun sensor. Extensive 
sun shielding and baffling, combined with the relatively 
small optical apertures of the Spider telescopes, allows 
the instrument to scan close to the Sun for increased 
sky coverage. 

The design of Spider has been extensively optimized 
to take full advantage of the low millimeter-wave back- 
grounds available from a stratospheric balloon platform, 
as well as to ameliorate polarized systematics to the level 
needed to characterize .B-mode polarization. Spider em- 
ploys detector arrays very similar to those successfully 
fielded in Bicep2 and the Keck array, but tuned for much 
lower noise- equivalent temperatures. Extensive filtering 
and cold (< 3 K) baffling within each instrument greatly 
reduces stray photon loading on the bolometers. Spi- 



der's simple, telecentric optics limit the contribution of 
polarized optical sidelobes. A 4 K half-wave p late at the 
aperture of each telescope ( Bryan et al.|20 10) is stepped 
in orientation periodically throughout the observation 
period in order to modulate the polarization signal of the 
sky with respect to any polarized instrumental systemat- 
ics. These systematics have already been characterized 
extensively during optical testing of the first Spider tele- 
scopes and through operation of similar technology at 
the South Pole. Tolerances and achi eved performance 
are described in more detail in Section l2~2l 



We refer t he re ader to Filippin i et al. l ( 2010 ) , Runyan 
(2010), and Gudmimdsson et al. l |2010[ ) for a more 



et al 
detai 



ed description of the Spider instrument 
2.2. Systematics Review 



Precise control of instrument systematics is crucial for 
achieving Spider's science goals. Table [2] summarizes 
our current understanding of a variety of sources of sys- 
tematic error, including gain uncertainty, pointing and 
beam effects, half-wave plate nonidealities, and sensitiv- 
ity to the Earth's magnetic field. We characterize each 
systematic effect using a suite of simulations with no in- 
put -B-mode power, by assuming a target level of control 
over the relevant parameters, and measuring the level of 
the resulting false -B-mode signal at t = 100. The design 
requirement is to reduce systematic error to a level of 
(43 nK) 2 , the value of the primordial r = 0.03 B-mode 
power spectrum at I = 100 in units of £(£ + l)Ce/ (2n). 

In most cases, the target values for each parameter are 
derived from a wel l-established simu lation pipeli ne de- 
scribed in d etail in |MacTavish et al!| ( |2008[ ) and O'Dea 
et al. (20111. The simulations include a detailed model 
ot polarized Galactic dust emission, which interacts non- 
trivially with instrumental effects. The aim of these sim- 
ulations is to quantify the extent to which a false B- 
modc signal is produced from systematics that induce 
/ — > Q, U or Q f-> U mixing, when no attempt is made 
to correct for them. These simulations assume a scan 
strategy and observing region that differ from the cur- 
rent baseline. However, internal simulations done using 
the McMurdo observing strategy described in Section [4~l] 
suggest that science requirements will be comparable. 
Further simulations are ongoing. 
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TABLE 1 

Spider observing bands, pixel and detector counts, and single-detector and single- telescope FPU sensitivities 



Band Center 
(GHz) 


Bandwidth 
(GHz) 


Beam FWHM 
(arcmin) 


Number of 
Spatial Pixels 


Number of 
Detectors per FPU 


Detector Sensitivity 

( /LtKoMBV^ ) 


FPU Sensitivity 
( /uKombVs ) 


90 


22 


49 


144 


288 


150 


10 


150 


36 


30 


256 


512 


150 


7 


280 


67 


17 


256 


512 


380 


18 



Note. — Each FPU sensitivity is obtained by dividing the corresponding single-detector sensitivity by V ^dct > assuming 
a detector yield of 85%, slightly below the average of the delivered focal planes. The total experimental map depth at each 
frequency scales inversely as the square-root of the number of FPU ■ flights for that frequency. The quoted sensitivities at 
90 GHz and 150 GHz are our current best estimate based on in-situ measurements of signal and noise using an aperture filling 
4 K load. The 280 GHz sensitivity is scaled from the average in-flight sensitivity of BOOMERANG at 245 GHz and 345 GHz. 



TABLE 2 

Expected sources of systematic error for Spider 



Source 


Target 


^residual InBB 

Woo /'-100 


Current Status 


Relative gain uncertainty: A(gi / 52) / (<7i / 32) 


0.5% 


17% 


0.1% Achieved by BOOMERANG 3, 


Differential pointing: (ri — r2)/u 


5% 




1% Measured in Spider 


Differential beam size: (a\ — CT2)/<7 


0.5% 


| 20%+>* 


0.3% Measured in Spider 


Differential ellipticity: (ei — E2)/2 


0.6% 




0.15% Measured in Bicep2 


Absolute polar angle calibration: A'J'abs 


1° 


17% 


0.7° Achieved by BlCEP b 


Relative polar angle knowledge: A*I/ rc i 


1° 


6% 


0.1° Achieved by BlCEP b 


Telescope pointing uncertainty: Ab 


10' 


6% 


2.4' Achieved by BOOMERANG 


Beam centroid uncertainty: Ac 


1.2' 


12% 


Achieved by BlCEP b 


Polarized sidelobes (150 GHz): G s 


-17 dBiH 


8%* 


Achieved by BlCEP b 


Optical ghosting: G T /Go 


-17 dB 


6% 


Achieved by BlCEP2 d 


HWP differential transmission: p 


0.7% 


10%t 


Achieved by Spider c 


Magnetic shielding at focal plane: <& s 


10 ^K CUB /B e 


3%t 


Achieved by Spider £ 



Note. — The target for each parameter is set so that the maximum false B-mode signal is less than the B- 
mode power spectrum at £ = 100 for r = 0.03. The third column provides the expec ted level of residual fo r each 
systematic effect if the target is exactly met (unless otherwise noted, numbers are from |MacTavish et al.|2008|) . The 
conventions for beam-related systematic errors foll ow those in [Bock et al.| |2008[ l. B e is tne strength ot tne iiarth's 
magnetic field. The optical ghosting simulations of MacTavish et al. ( 2008 ) primarily probe E — > B contamination 
and not the potentially more important T —> B contribution that would be expected given nontrivial polarization 
of the optical ghosts. We are currently revisiting the ghosting simulations to incorporate our know ledge of their 
polarized response. Given the frequency dependence of the foreground emission (see Section |3.1| l, requirements 
for polarized sidelobes at 150 GHz are sufficient for observations at 90 and 280 GHz. The Current Status column 
indicates whether each target is met given current measurements (for Spider, this corresponds to lab studies). 
When no reference is given, the current status is based on unpublished measurements. Adding all systematic errors 
in quadrature gives a total false signal of ~37% of the r = 0.03 B-mode signal at 

II Gain at 30° from t he telescope boresight 
t |Q'Dea et"al] ( |201l| ; * This paper (see Section 
Amplitude of the resi dual for the three sourc es of systematic error c om bined 

I Jones et al. | ( |2006a[ | ; 



100. 



Takahashi et al.| (|2010|); 



|Runyan et al.| | |2010[ | 



d |Aikin et al. 



(2010) 




Fig. 2. — Left: Polarized sidelobe profiles measured at 150 GHz in BlCEP (green) and models used in the analysis presented in Section |2,3| 
We consider two models: a sidelobe profile comparable to Bicep's (blue), and Spider's target profile, which assumes an additional 10aB 
attenuation starting at 12° from the telescope boresight (red). Right: B-mode power spectra of the sidelobe pickup obtained with the 
sidelobe profiles shown in the left panel. The r = 0.03 primordial B-mode power spectrum is shown for comparison. The contamination 
is well below the cosmological signal at i = 100 for both profiles. However, the BlCEP-like model leads to substantial contamination at 
large angular scales (£ < 25). Spider's target profile ensures that the false B-mode signal is comparable to or fainter than the cosmological 
signal at low multipoles, and negligible at I ~ 80, where the primordial B-mode power spectrum peaks. 
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Many of the targets listed in Table [2] for systematics 
related to design elements common to Doth Bicep and 
Spider have been met by the former ( |Takahashi et al. 
20101, which ensures that Spider is in a position to do 
so as well. Additional sources of systematic error due to 
half-wave plate nonidealiti es and magnetic fie ld pickup 
are investigated in detail in O'Dea et al. (2011 ). In both 
cases, we find that the level of if- mode contamination is 
small compared to the r = 0.03 inflationary signal, even 
without special effort to correct for these systematics. 

Spider's control of systematic errors is further en- 
hanced by the use, with each telescope, of a rotating 
half-wave plate that is periodically stepped so that the 
same patch of sky is observed at multiple orientation 
angles, thereby augmenting the natural angular modu- 
lation due to the motion of the sky. This homogenizes 
the statistical noise in polarization and provides pow- 
erful rejection of a variety of instrumental systematics. 
By locating the half-wave plate skyward of the telescope 
aperture, we ensure that only the sky signal is modu- 
lated and that the effects of beam mismatches between 
detector pairs in each spatial pixel are greatly reduced. 

Gain drifts contribute significantly to the systematics 
budget. Boomerang achieved a relative gain uncer- 
tainty less than 0.1%, a level fi ve times smaller than th e 



target established for Spider in MacTavish et al. (2008), 



using an on-board calibration lamp (|Masi et al 



2006) 



We expect Spider to reach a similar level of control us- 
ing regular modulations of the detector bias currents. 
This technique is currently undergoing laboratory tests. 
A BoOMERANG-style calibration lamp is also an option, 
pending the results of these studies. 

Thermal fluctuations within the instrument will in- 
evitably create a signal in the Spider bands. Changes 
in focal plane temperature (e.g., with payload altitude) 
can mimic a sky signal by changing the power flow 
through the bolometer support membranes. Spider's 
TES bolometers self-heat to a temperature largely inde- 
pendent of that of the focal plane, rendering them less 
sensitive to such fluctuations than a comparable semi- 
conductor bolometer. Since the monolithic fabrication 
of each detector array ensures good uniformity in ther- 
mal response within each polarization pair, such fluc- 
tuations do not propagate as strongly into a false po- 
larization signal. A passive thermal filter, consisting of 
stainless steel heat capacity blocks between the detectors 
and thermal strap, supresses any thermal fluctuations on 
time scales comparable to the scan (by a factor of ~ 20 dB 
at 30 mHz). 

Since about 25% of each beam terminates on the op- 
tical stop, changes in the stop temperature will also 
mimic a sky signal. Essentially all such internal side- 
lobes terminate on the cooled (~ 1.8 K) optics sleeve, 
however, whose low temperature limits the sensitivity of 
the bolometers to such drifts. 

We report on simulations specifically designed to eval- 
uate the level of the contamination expected from polar- 
ized sidelobes in the next section. 

2.3. Polarized Sidelobes 

Coupling through the beam sidelobes to emission from 
the Sun, the Earth, and the Galactic plane will produce 
a spurious polarized signal. In order to establish mission 
requirements for Spider's sidelobe response, we convolve 



a polarized sidelobe profile with a temperature map of 
the Galactic emission at 15 GHz based on the sky model 
developed in |Q'Dea et aL (20111 and discussed further 
in Section |3.2[ Because of its complex morphology and 
proximity to spider's observing region, the Galaxy is in- 
deed expected to represent the most significant source of 
contamination on the scales of interest. Spider's optical 
system is a cryogenic refracting telescope, based on that 
developed for Bicep. Unlike Bicep, Spider employs a 
conical baffle with its base located inside the cryostat 
and extending roughly 0.5 m out from the surface of the 
vacuum vessel. The rim of the baffle is asymmetric to 
prevent direct illumination of the half-wave plate by the 
balloon and the Earth's limb, and intercepts the field 
of view between 10°-12°, in contrast with the 15° cut- 
off angle for the Bicep telescope forebaffie. Spider's 
main beam and near (2°-12°) sidelobe profiles are in- 
formed by a physical-optics model of the ideal optical 
chain internal to the cryostat, which includes both lenses 
and the half-wave plate. Accurately modeling Spider's 
far (> 12°) sidelobes requires a complete simulation of 
the optical system including sun shield, baffle, gondola, 
and balloon that has yet to be performed. Instead, we 
first model Spider's far sidelobe profile as a power law 
out to 50°, where the beam is truncated. This model is 
consistent wit h measurements of Bic ep 's polarized side- 
lobe response ( Takahashi et al.|20 l6|, which is represen- 
tative of Bicep s relatively less restrictive forebaffie. The 
profile we use is shown in the left panel of Figure [2j 

The £?-mode power spectrum derived from the sidelobe 
model described above is shown in the right panel of 
Figure [2j Compared to the primordial r — 0.03 -B-mode 
signal, it is over 12 times fainter at I = 100, while it is 
significantly brighter at large angular scales (£ < 25). 

Unsurprisingly, the Bicep forcbafflc, which was de- 
signed to provide rejection at the level of r ~ 0.1 while 
observing at relatively high Galactic latitudes, does not 
provide sufficient rejection of the foreground signal for 
Spider's large-scale measurements. Our simulations 
suggest that Spider requires an additional ~ 10 dB of 
attenuation beyond 12° off axis in order to reduce the sys- 
tematic contribution to a level near or below the r = 0.03 
-B-mode spectrum at all scales of interest. The contam- 
ination is then over an order of magnitude fainter than 
the cosmological signal at I ~ 30, and entirely negligi- 
ble at the I ~ 80 peak (see Figure [2]). Although further 
measurements are needed to demonstrate that Spider's 
baffle will provide the required factor of ten improvement 
over the Bicep measurement, it appears to be well within 
reach. For instance, geometrical theory of diffraction cal- 
culations indicate that the introduction of an additional 
low angle of incidence diffraction edge to the baseline baf- 
fle provides the extra —10 dB rejection req uired beyond 



a few degrees of t he geometrical shadow (|Keller||1962 
Born fc WolffT999l ). 



3. GALACTIC FOREGROUNDS 

3.1. Polarized Foregrounds in the Spider Field 

During its first two flights, Spider will map the mi- 
crowave sky at frequencies ranging from 90 GHz to 
280 GHz. At 90 GHz, the large-scale {£ < 10) polarized 
emission from interstellar dust is expected to be at least 
an order of magnitude brighter than the r = 0.03 primor- 
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Fig. 3. — Left: A multipole-by-multipole comparison of the levels of statistical noise at each frequency in Table Rafter both flights and 
of the astrophysical and cosmological signals, including Galactic foregrounds, lensing, and the B-mode signal assuming a scalar-to-tensor 
ratio r = 0.03. Right: Level of Galactic foreground emission in the nominal Spider field (/ s k y = 10%) compared to constituent trial fields 
with / s k y = 2%. The optimal 2% of the sky, outlined in Figure Hi has the lowest polarized dust emission at 50 < £ < 200 not only within 
the Spider region, but also within the entire area of sky accessible from a McMurdo LDB flight. Spectra are also shown for a trial field 
centered on (a = 0°, S = —57°), which lies within the "Southern Hole," a target region used by several ground-based experiments. Note 
that our foreground model does not include po lariz ed synchrotron emission from the Galaxy given the weakness of this signal in our bands 
in the Spider region (see discussion in Section|3.1[l. 



dial B-mode si 



gnal that Spider aims to detect (Dunk- 



ley et al. 2009). At smaller scales, the angular power 



spectrum ot the polarized dust emission appears to be 
compatible with the pow er-law dependence C^ ust oc £~ 2 - 6 
(Cho & Lazarian 2010), which results in its amplitude 
being comparable to that of our target B-mode signal at 
£ ~ 40. Models indicate that this multipole dependence 
might break down in favor of a steeper decline at mul- 
tipoles higher than ~ 10 3 . Assuming that the polarized 
emission from interstellar dust obeys the same frequency 
dependence as its total emission, these results can be ex- 
trapola ted to higher frequencie s using the best-fit FDS 
model (Finkbeiner et al. 1999). The resulting angular 
power spectra are ~ 9 and ~ 5 x 10 2 times higher than 
the 90 GHz signal at 150 GHz and 280 GHz, respec- 
tively. Throughout this paper, we assume spinning dust 
emission to be unpola rized, consistent with theore tical 
expectations (see, e.g., Lazarian & Finkbeiner 2003). 

Among other Galactic foregrounds, only synchrotron 
emission is expected to be significantly polarized at mi- 
crowave frequencies. We use the WMAP 23 GHz data 
to evaluate the amplitude of the polarized Galactic syn- 
chrotron emission in the Spider region at frequencies 
ranging from 90 GHz to 280 GHz. The extrapolated 
signal is negligible both at 150 GHz and at 280 GHz 
compared to the r = 0.03 B-mode signal. At 90 GHz, 
its large-scale power spectrum is a factor of two higher 
than that of the target signal, but a factor of five lower 
than that of the polarized emission from interstellar dust. 
We find the multipole dependence of the polarized syn- 
chrotron emission to be well described by the power law 
2 5 (compa tible with WMAP's full-sky 
Page et al.|2007 ), which brings this spectrum 
ot the r = 0.03 B-mode signal by I ~ 30. At 



^synchrotron £ 



estimate in 
under that 

the I ~ 80 B-mode peak, this foreground is an order of 
magnitude fainter than the target signal. As a result, we 
do not expect the Galactic synchrotron emission to limit 
Spider's ability to detect B-modes at the r — 0.03 level. 
Since the Galactic plane will be masked during the cos- 



mological analysis, we do not consider the polarization of 
the free-free emission that might be in duced by Thomson 
scattering at the edges of Hn clouds ( Rybicki & Light' 
. Free-free emission is otherwise unpolarized. 
le considerations in this section, we expect the 



man||1979 
Given t 



polarized emission from interstellar dust to be the one 
Galactic foreground that will hinder Spider's ability to 
dete ct a primordial B-mode signal. We show in Sec- 
tion [3]2] that it is possible to select observing fields acces- 
sible irom a McMurdo flight in which the polarized emis- 
sion from interstellar dust is up to an order of magnitude 
fainter than discussed above at all Spider frequencies. 
However, even with this reduced level of contamination, 
Galactic dust is still the dominant polarized sky signal 
at and above 90 GHz. Section [5] therefore addresses the 
question of how Galactic emission and cosmological sig- 
nal can be separated given an appropriate combination 
of frequency bands. We describe the Galactic foreground 



model required for both these analyses in Section 3.2 



3.2. Model for the Polarized Emission from 
Interstellar Dust 

A model of the sky that Spider w ill observe was pre 



sented in detail in O'Dea et al. (2011 ). In this section, we 



summarize the procedure followed to create a template 
of the polarized emission from interstellar dust that takes 
into account the three-dimensional Galactic structure. 

The nonsphericity of interstellar dust grains and their 
ability to align with the Galactic magnetic field (here- 
after, the Field) were first put forth in an effort to 
explain t he polarization of starli ght observe d indepen- 



Hall 



(jl949j). With its 
icular to the Field 



dently by Hiltner (19491 and by 
long axis prcfcrcntially "aligned perp enc 
(see, e.g., |Hoang &: Lazarian|[2008| , a dust grain absorbs 
more radiation in the direction perpendicular to the local 
Field line than in the direction parallel to it. This re- 
sults in a net linear polarization of the incident starlight 
in the direction parallel to the Field. The polarization 
of the emission from dust grains is the counterpart of 
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TABLE 3 

Summary of general information on the Spider mission 



Launch Location 


McMurdo Station, Antarctica 


Launch Date 


12/2012 (Flight 1), 12/2014 (Flight 2) 


Flight Duration (target) 


20 days per flight 


Altitude (target) 


36,000 m 


Flight Path 


Circumpolar, typically 73° S < latitude < 82° S 


Sky Coverage 


/ sky ~ 0.1, 10 < I < 300 



Note. — The flight schedule provided in this table is consistent with the 
state of hardware development as of June 2011. 



this absorption phenomenon. As a result, the emission 
from interstellar dust is expected to be polarized per- 
pendicularly to the sky-projected direction of the Field. 
Archeops and WMAP both detected the polarization of 
the Galactic dust emission at microwave frequencies with 
high significance (Benoit et al. 2004 Page et al. 2007), 
and found degrees of polarizatio n over the sky compatibl e 
with theoretical expectations (Draine & Fraisse 2009). 
However, neither experiment had the combination of res- 
olution and sensitivity necessary to produce a map of this 
emission usable by Spider. As a result, we must rely on 
modeling to estimate the characteristics of the polarized 
emission from interstellar dust in the Spider held. 

Despite tremendous recent progress in our understand- 
ing of the alignment of interstellar dust with the Gal actic 
magnetic field (see, e.g., the review in |Draine||2011[ ), we 
are still far from a full theoretical understanding or this 
physical process. Modeling the polarized emission from 
Galactic dust therefore requires simplifying assumptions. 
It is customary to assume (i) that each grain is aligned 
with the Field with its long axis exactly (as opposed to 
preferentially) perpendicular to the field line; (ii) that, 
given a field strength, each grain has the same polariz- 
ing effect; and (iii) that the degree of polarization of the 
dust emission induced by a grain is proportional to the 
square of the magnetic field strength at the location of 
the grain. With these assumptions, and given a three- 
dimensional model of the dust distribution and of the 
Galactic magnetic field, we compute a full-sky template 
of the polarization fraction and angle of the Galactic dust 
emission. Details of the three-dimensional models and 
of the calc ulation performed t o derive the template are 

Since the magnetic field 
urbulent component, the 



provided in |Q'Dea et al.| ( |2011 [) 
model includes a small-scale 
resulting depolarization effect is taken into account. The 
overall amplitude of the template for the polarized frac- 
tion of the Galactic dust emission is a free parameter, 
which we set to match the average value of 3. 6% de 



rived by W MAP outside of the Galactic plane (Kogut 
et al.|2007 ) . We obtain polarized intensity maps by mul 



tiplyi ng~the template b y the dust intensity map from 
Finkbeiner et al. (TT999), using their model 8 to account 
for its frequency dependence . 

The left panel of Figure [3] shows the power spectra of 
the polarized Galactic dust emission in the Spider field 
for all three frequencies in Table 111 along with the per- 
multipole statistical error for each band after both flights 
as predicted by a Fisher analysis (see Section [5]). At a 
given frequency, the brightness of the polarized emission 
from interstellar dust in the Spider field is compara- 
ble to that of the full-sky average of this signal (Galac- 
tic plane excluded) , whose characteristics were described 



in Section |3.1| Interestingly, it is also comparable to 
that of the polarized Galactic dust emission in a typical 
/sky = 2% patch in the popular "Southern Hole" region, 
as shown in the right panel of Figure[3j even though Spi- 
der will observe five times as much sky as covered by this 
patch. This appears to indicate that although the South- 
ern Hole is believed to be the region of the southern sky 
most free of dust emission, it might not be the region 
most free of polarized dust emission. Finally, it is worth 
noting that the Spider region encompasses the cleanest 
2% of the sky accessible from a McMurdo LDB flight, 
and that a large majority of its component fields exhibit 
significantly less polarized Galactic dust emission than 
the region average (the relevant power spectra are also 
shown in the right panel of Figure [3]) . This will provide 
valuable cross-checks to evaluate the level of foreground 
contamination in the Spider maps. 

4. OBSERVING STRATEGY 
4.1. Overview 

Spider will launch its first flight from McMurdo Sta- 
tion in Antarctica in December 2012, and observe the 
southern sky for about 20 days. The flight parameters 
are summarized in Table [3j The observing strategy is 
designed to maximize sky coverage and crosslinking over 
the cleanest regions of the sky accessible from a Mc- 
Murdo LDB flight. Figure [4 shows the intensity and 
polarization of the Galactic g ust emission in the Spi- 
der region ( O'Dea et al.||20lT |, as well as the relatively 
foreground-free CMB sky in the 94 GHz WMAP band. 

The baseline scan strategy consists of azimuthal scans 
punctuated by elevation steps and rotations of the half- 
wave plate. The scan velocity profile is sinusoidal with a 
maximum acceleration of 0.8 degrees/s 2 (for a scan pe- 
riod of ~ 45 seconds). The exact scan amplitude and 
center is adapted at each elevation step to (i) avoid ob- 
serving within 90 degrees of the Sun, as set by the design 
of the sun shield, and (ii) maintain the instrument bore- 
sight pointing at the region of the southern sky outside 
of the Galactic plane (0 to 80 degrees right ascension). 
Elevation steps occur every hour, in sync with local side- 
real time (LST). We reach maximum elevation (40 de- 
grees) at 12:00 LST, and minimum elevation (28 degrees) 
at 24:00 LST. By syncing the elevation steps with side- 
real time, we ensure maximal sky coverage in declina- 
tion. The sky coverage repeats each sidereal day, provid- 
ing both redundancy and a rich set of consistency tests. 
Once per day, when sky rotation is at a minimum at local 
midnight or noon (depending on when observations be- 
gin) , we rotate the half- wave plate of each instrument by 
22.5°, thereby switching Q to U in the instrument frame. 

We set aside time after each half-wave plate rotation 
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Fig. 4. — Spider observing region and scan profile. All four panels show the same portion of the sky, in equatorial coordinates, smoothed 
with a 30' beam. The southern Galactic pole (black +) is overplottcd, along with the 10- and 20-degree Galactic latitude lines (dashed). 
Top left: IRAS 100 /im data showing dust morphology, the Spider observing region (white outline), and the south ecliptic pole (white *). 
Also shown are the BOOMERANG and BlCEP fields (left and right gray outlines, respectively), and the region of minimum foreground 
contam inat ion in the Spider field (black outline). Bottom left: Polarized dust emission amplitude at 150 GHz, according to the model in 
Section [3.2| Clear differences in dust emission morphology are visible between this model and the IRAS data. Top right: WMAP 94 GHz TT 
data in the same area, showing relative absence of foreground contamination in the Spider observation region. Botto m r ight: Distribution 
of integration time, averaged over all detectors in a single 150 GHz focal plane, for the observing strategy in Section |4.1| 
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Fig. 5. — The Spider scan profile. All three panels show the same portion of the southern sky, in equatorial coordinates, smoothed with 
a 30' beam. Top: A 24-hour period for a single detector, illustrating the change in telescope orientation throughout the day. Five-minute 
periods every three hours are shown. Detector pairs in a single p ixel are oriented 22.5 and 112.5 degrees relative to the scan direction. 
Bottom left: Fraction of excess variance (see definition in Section |4.1| due to anisotropic angular coverage in the Stokes Q and U maps 
for the above detector over a four-day period. Bottom right: Fraction of excess variance for a full focal plane of 512 bolometers, over a 
four-day period. This observing profile covers 9.5% of the sky, of which 96% is observed with near-isotropic coverage in crossing angles to 
reconstruct the Stokes Q and U parameters with < 1% excess variance. 
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to perform beam and pointing calibrations on Galac- 
tic sources. Every half hour, near the beginning and 
end of each constant-elevation period, we perform a brief 
(~ 5 s) in-flight measurement of detector responsivity by 
slightly stepping the detector biases. Each sub-Kelvin 
cooler must be re-cycled every 72 hours. To reduce the 
instantaneous load on the 1.6 K stage, these cycles must 
be staggered. The cycling of one instrument has no im- 
pact on the others, allowing for a ~ 85% observational 
duty cycle. The process of cycling all six coolers requires 
four hours for all instruments to return to the nominal 
base temperature. 

We estimate the effect that the observing profile will 
have on the Stokes T/Q/U reconstruction, assuming an 
ideal case of no beam asymmetry, no T — > Q, U leakage 
(i.e. perfect polarization efficiency) , and no noise correla- 
tion between detectors. We calculate the pointing matrix 
A p for each pixel p, as 



Ap — 



1 

cos 20 



cos 2-0 
cos 2 20 



sin 20 
cos 20 sin 2^ 



sin 20 cos 20 sin 20 sin 20 



(1) 



where (. . .) indicates the average over all i = 1 . . . N b S:P 
observations along p at the series of polarization angles 
{4>p,i}- The inverse of A p gives an estimate of the signal 
T/Q/U covariance for each pixel. In the ideal case, we 
would have A~ x = diag(l,2,2). We therefore construct 
a figure of merit, which we call the "fractional excess 
variance in polarization," from the QQ and UU elements 
of Al 1 : 



1 



'QQ 



-V2 



>uu 



V2 



-i 2 



(2) 



This figure of merit is shown in Figure [5j for a single de- 
tector, and for a full focal plane of 256 pixels (2 detectors 
per pixel). The excess variance is below 1% over 96% of 
the observed region. We will fly with pairs of telescopes 
clocked at 45° relative to each other, which allows a si- 
multaneous measurement of both Q and U parameters, 
significantly improving crosslinking in the maps and fur- 
ther reducing the excess variance. 

4.2. Filter Transfer Function Fn 

The interplay between scan strategy and instrumental 
and environmental effects plays a significant role in all 
CMB experiments. These instrumental effects include 
scan-synchronous artifacts, ground-fixed signals, atmo- 
spheric contamination, and low-frequency instability in 
the instrument. Signal estimation pipelines tend to be 
tailored to each experiment (or flight of the same exper- 
iment) according to the particular combination of effects 
that are most relevant to each dataset. 

While these algorithms, which include template re- 
moval, common mode decorrelation, and (statistical) 
least-squares mapmaking, are generally quite effective at 
removing systematic effects from the data, they never- 
theless result in a degradation of sensitivity and the loss 
of fidelity to particular spatial modes of the astrophysical 
signal relative to an idealized dataset consisting only of 



the signal and Gaussian uncorrelated white noise. Spi- 
der is no exception, and we study the impact of these ef- 
fects through time domain simulations of the experiment 
that capture as fully as possible the range of instrumental 
and environmental effects that are present in the data. 

The impact on the fidelity of the recovered CMB power 
spectrum can be approximately characterized via a cal- 
culation of the transfer function Fg 1 as determined by the 
comparison of a modest ensemble of signal simulations to 
the input sky ( Hivon et al. 2002 1 . An understanding of 



this transfer function is necessary to relate NET esti- 
mates of any experiment, for example those in Table [T] 
to errors on CMB power spectra or experimental upper 
limits. As we will show, the relatively benign strato- 
spheric balloon environment coupled with Spider's large 
sky coverage provide polarized maps of the sky with ex- 
tremely low noise and high fidelity to angular scales rang- 
ing from the scale of the map to that of the beam. 

For Spider, we expect the large angular scales to be 
dominated by noise at time scales below the half-period 
of an azimuthal scan, which varies between ~ 20-25 s, 
depending on the scan amplitude. Contributions to the 
noise include a scan synchronous component from the 
residual atmospheric column as well as 1// noise from 
instability on the bias and instrumental backgrounds. 
While the stationary Gaussian component of the noise 
can be optimally accounted for by the use of a least- 



squares mapmaker (see, e.g., |Jones et al.||2007| , we ex- 
pect, based on our experience with past balloon exper- 
iments, the need to apply additional filtering to remove 
a scan-synchronous component and reduce the impact of 
uncertainty in the low-frequency transfer function. 

We measure the net effect of signal processing on the 
power spectrum by calculating the transfer function F i 



following the formalism developed in|Hivon et al. (20021 



We calculate the power spectrum Cf 11 ' 1 of a signal-only 
CMB realization that has been smoothed by the beam 
(30' FWHM for the 150 GHz instrument). We then 
apply the observing profile shown in Figures [4] and [5] 
to generate noise-free time-ordered data for processing. 
These simulated data are then input to the signal esti- 
mation pipeline, which incorporates our estimate of the 
noise spectra, flagging and any additional filtering. We 

then calculate the power spectrum Cg^ of the resultant 
maps. Both the input and output maps are weighted by 
the distribution of integration time shown in Figure [4] to 
ensure identical cut-sky treatment. The transfer function 
is then simply the ratio of the input and output power 
spectra, averaged over an ensemble of realizations. 

Figure [6] shows the published transfer functions for Bi- 
CEP and Boomerang03, as well as two estimates for 
Spider that correspond to optimistic and pessimistic 
levels of scan-synchronous noise. Results are shown for 
the temperature anisotropy power spectrum only; mea- 
surements of the polarized filter function are ongoing, but 
are not qualitatively different. We multiply the transfer 
function with the beam to illustrate the combined ef- 
fect on the power spectrum. At large angular scales, we 
lose modes due to interplay between the filtering and 
the scan, while small angular scales are downweighted 
by the beam. The relatively severe required filtering ap- 
plied in Bicep reduces the product of the transfer func- 
tion and beam window function to nearly 50% at the 
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peak. In the optimistic scenario for Spider, we set the 
high-pass filter cutoff just above the azimuthal scan fre- 
quency (25 mHz), which leads to the peak nearing 90% 
at t ~ 50 even with a naively binned map. More realisti- 
cally, we have set the cutoff near 50 mHz, corresponding 
to a half-period of the scan. This choice is informed by 
the data from Boomerang, and corresponds to the most 
stringent filtering that had to be applied to the polariza- 
tion data from the 2003 flight (towards the end of the 
flight, at the relatively low altitude of ^30 km). Spider 
preserves fidelity to larger scales (lower multipoles) than 
Boomerang simply because of the larger sky coverage 
and higher scan rate, and not because the Spider data 
are expected to be qualitatively more stable or well be- 
haved than the Boomerang data. While this more ag- 
gressive filtering reduces the transfer function at i ~ 10 
by over 50%, the peak at I ~ 70 remains relatively high 
at 85%. 
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Fig. 6. — An estimate of the filter transfer function Fj T and 
beam function Bi for Spider's 150 GHz instrument, show n along- 
side transf er functions for typical ground-bas ed (BlCEP. |Chiang| 
|et al.|20io) and balloon-borne (B00MERANGQ3, pones et al.]^UUbb| l 
experiments. The transfer function for Spider is estimated rrom 
an ensemble of signal realizations assuming time-domain filtering 
of signals below the azimuthal scan frequency. The black and 
blue curves were calculated from high-pass filtered signals naively 
binned into maps, while the red curve was obtained using an iter- 
ative least-squares mapmaker that took into account both a high- 
pass filter and the expected bolometer noise PSD. 

The covariance of the band power estimates derived 
from the data depends on the effect of the experiment on 
both sample and statistical variance. The sample error 
on the measured power spectrum Ci increases as 1/y/Tg 
as the effective number of observed modes is reduced. 
While this increase in the sample variance has no impact 
on the ability of an experiment to set upper limits on the 
-B-mode spectrum, the impact of the processing on the 
statistical variance does. 

To account for detector noise, we use an iterative least- 



squares mapmaker similar to that in Jones et al. (20071 
on the filtered data. The inverse-noise kernel is esti- 
mated from a typical noise power spectrum, measured 
in our system with no optical loading on the bolometers. 
The absence of (Gaussian and white) photon noise in this 
measurement provides an accurate estimate of the instru- 
mental noise. The significant features in the noise power 
spectrum are (1) a 1// knee at ~ 100 mHz, and (2) ex- 
cess noise at frequencies above ~ 10 Hz, likely due to 



the thermal architecture of the bolometer island, which 
reduces the asymptotic transfer function by ~5%. The 
iterative method improves the net transfer function by up 
to 30% relative to the naively binned case at the largest 
scales. Recent bolometer noise measurements under op- 
tical loading with improved bias stability show that the 
1// knee is likely to be closer to 25 mHz in flight, which 
should further improve the large-scale transfer function. 

5. BAND SELECTION & PROSPECTS FOR 
FOREGROUND REMOVAL 

In this section, we evaluate the ability of Spider to 
detect inflationary £?-modes as a function of the fre- 
quency distribution of FPU-flights (number of focal plane 
units multiplied by number of flights). In addition to the 
90 GHz, 150 GHz, and 280 GHz channels whose charac- 
teristics are given in Table [T] we consider the possibility 
of adding a 220 GHz frequency band. This frequency 
was chosen to be half way between Spider's 150 GHz 
and 280 GHz channels, thereby eliminating the frequency 
coverage gap that results from the frequency distribution 
in Tablep] The 220 GHz band would have the same num- 
ber of spatial pixels and the same number of detectors per 
FPU as the 280 GHz band. Its bandwidth and its beam's 
FWHM would be 53 GHz and 21', respectively. The sen- 
sitivity of each of its detectors would be 320 fiK CMBy /s, 
for an overall FPU sensitivity of 15 /iK CMBV /s. For pro- 
grammatic reasons, we require that observations at each 
frequency during each flight be performed by at least 
two FPUs. For each Spider flight, we therefore have the 
option to acquire data at up to three frequencies to be 
chosen among 90 GHz, 150 GHz, 220 GHz, and 280 GHz. 
Since, as discussed in Section[3j the emission from Galac- 
tic dust is expected to become stronger with increasing 
frequency, Galactic foregrounds are expected to be lower 
at 90 and 150 GHz than at 220 and 280 GHz. As a re- 
sult, we choose 90 GHz and 150 GHz as two of the three 
Spider bands for the first flight, and perform numerical 
simulations to distribute the remaining two FPUs. 

We simulate a hundred maps at each frequency with 
independent CMB and noise realizations. CMB maps 
are obtained with WMAP's best fit ACDM cosmolog- 
ical parameters, r = 0.03, and CMB lensing included. 
We approximate the noise in Spider as white, governed 
by the parameters in Table [l] (for the 220 GHz band, 
we use the parameters listed m the previous paragraph), 
and weighted by the expected number of o bser vations per 
pixel for the observing strategy in Section |4~Tj Each sim- 
ulated map is added to a template for the Galactic dust 
emission at the appropriat e wavelength Dust t emplates 
are extrapolated from the O'Dea et al. (2011) 90 GHz 
map described in Section [3 assuming a single power-law 
dependence on frequency with spectral index /3<j = 1.7, 
the latter bei ng the average o f the o ne-index best fit 
to model 8 of |Finkbeiner et al l (jl999l and of W MAP's 
"base" fit dust spectral index (Gold et al.||2009|). Note 



that we do not attempt to use a two-component dust 
model. In some of the cases we consider in this section, 
this would indeed lead to more parameters than we can 
possibly fit given the number of observing bands in play. 

With two pairs of FPUs occupied by either 90 GHz or 
150 GHz detectors, we use a simple pixel-by-pixel least- 
square foreground fitting procedure to evaluate Spider's 
first flight's ability to detect B- modes when the two re- 
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TABLE 4 

Spider FPU frequency distribution and per-band cumulative noise 



Flight 


FPU Distribution 


Cumulative Noise (/iK C MB/deg 2 ) 
90 GHz 150 GHz 280 GHz 


Spider 1 
Spider 2 


3 X 90 GHz; 3 X 150 GHz 
2 X 90 GHz; 2 X 150 GHz; 2 X 280 GHz 


0.27 0.20 

0.21 0.16 0.62 



Note. — The multiple telescopes at each frequency are rotated with respect to one 
another to simultaneously recover Q and U in each daily map. The distribution of 
observing frequencies is chosen following the procedure described in Section [5] so as to 
maximize Spider's ability to detect inflationary B-modes at angular scales corresponding 
to t ~ 80. The cumulative noise numbers use the parameters listed in Table [T] 



maining FPUs are assigned to (i) 90 GHz and 150 GHz, 
(ii) 220 GHz, or (hi) 280 GHz. In order to have similar 
signal-to-noise ratios at the two low- foreground frequen- 
cies, we do not consider cases in which the numbers of 
FPUs at 90 and 150 GHz differ. We model the data as 



S — S CMB + Sd,^ 



(3) 



where S is the usual set of Stokes parameters, the index 
d refers to dust, and Vq = 90 GHz, and fit the seven 
free parameters to the simulated data. In all cases, we 
include as part of the data a Planck 217 GHz map sim- 
ulated in the same way as the Spider maps with the 
instrum ent characteristics published in the Planc k early 



papers (Planck HFI Core Team: Ade et aL|2011|). 

Given Spider's focus on characterizing the i ~ 80 B- 
mode peak, we use the £ = 80 .B-mode signal-to-noise 
ratio (SNR) computed from residual maps as a figure of 
merit to assess the quality of the CMB reconstruction. 
Among the three cases we consider for the first Spider 
flight, having three FPUs at 90 GHz and three FPUs at 
150 GHz leads to the highest B-mode SNR at £ = 80, 
and therefore to the best reconstructed CMB map for 
our purposes. This configuration turns out to be also 
favored from a detector development point of view. 

After the first Spider flight, we expect the per- 
multipole statistical error at 150 GHz to be at a level 
comparable to that of the power spectrum of our fore- 
ground model at t ~ 80. With Galactic dust now limit- 
ing our ability to detect inflationary i?-modes, we seek to 
increase Spider's frequency coverage in order to gather 
multi-frequency foreground information, which will help 
constrain the model in Equation ([3]). Given the pro- 
grammatic constraints previously mentioned and con- 
tinued observations at 90 GHz and 150 GHz, this can 
be achieved only by flying a pair of either 220 GHz or 
280 GHz FPUs. Performing an analysis similar to that 
described above for the first Spider flight, we select the 
280 GHz band over the 220 GHz channel, which leads 
to the FPU • flight distribution listed in Table 13] along 
with the cumulative noise in each band. Figure R] shows 
the resulting per-multipole statistical error in each band 
after both flights, as well as the overall statistical error 
for the full mission, both derived from a Fisher analysis. 
The expected Planck HFI statistical error from a similar 
analysis is also shown for comparison. 

Figure [8] shows forecasts for the constraints on the 
tensor-to-scalar ratio r from Spider (after one and both 
flights) and from Planck. These constraints are de- 
rived using the method described in detail in Farhang 
et al. (submitted), which we summarize below. 




10 Multipole / 



Fig. 7. — Theoretical E-mode and S-mode (EE and BB) 
power spectra and projected per-multipole statistical errors for 
Spider after both fl ights and tor Planc k HFI. E-mode detec- 
tions from W MAP (Larson et al. | |2011), BlCEP JChiang et al.| 
|2010|), QUIET {Q uiet Collaboration: ~Bischoff et aT.||ijUl(J[ |, ancT 
QUaL) ( |Brown et al.|[2009[ l are also shown tor comparison. The 
B-mode spectrum is the sum of two components (dotted lines): in- 
flationary gravitational waves, shown here for r = 0.03, and weak 
gravitational lensing of B-mode polarization. The noise curves for 
Spider and Planck are derived from a simple Fisher analysis, as- 
suming no foreground contamination. It is likely that foreground 
emission will limit Planck's measurement of the BB bump below 
t ~ 10 (corresponding to an angular scale of > 20 degrees) induced 
by reionization. Spider, is optimized to cover the I ~ 80 peak in 
the gravitational wave spectrum. The Spider band powers have 
been truncated at a multipole of seven as the limited sky coverage 
will prevent us from probing larger scales. 

Assuming a uniform a 'priori probability distribution 
for the cosmological parameters q, the a posteriori dis- 
tribution function's dependence on q is given by the like- 
lihood £a(<i) = -P(A|q), the probability of the data A 
given the parameters, such that 



21n£ A (q) 



A f C t -i A - 



In |(2tt) 



(4) 



where A is the full set of CMB temperature and polar- 
ization maps including both signal and noise, iV p i x is the 
number of pixels in the maps, and C to t = Cn + Qr(q) 
is the theoretical pixel-pixel covariance matrix with con- 
tributions from both the parameter-dependent signal co- 
variance Cr(q) and the generalized noise covariance Cn- 
The latter includes uncertainties from foreground sub- 
traction as well as from noise in the maps. 

Here, we calculate the likelihood in Equation ([JJ on a 
two-dimensional grid consisting of the cosmological pa- 
rameters r and r, the latter referring to the Thomson 
scattering optical depth to reionization. Although extra 
dimensions could be added, Farhang et al. (submitted) 
show that r is only weakly correlated with other stan- 
dard cosmological parameters in the small-r, small- / s k y 
limit of relevance to Spider. Since stringent constraints 
on t are an objective of large-scale experiments such 
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as Planck, we marginalize over it in order to evaluate 
whether insensitivity to it will impact Spider's ability 
to constrain r. It appears that with our sky coverage 
fixed at 10%, r and r are essentially uncorrelated. 

The results in Figure [8] assume / s k y = 0.75 for 
Planck and / s k y = 0.10 for Spider, and are given for 
two cases: "with" and "without" foregrounds. In the 
former case, foreground residuals after the removal pro- 
cedure described earlier in this section are included as 
contaminants, whereas in the latter case, it is assumed 
that there is no foreground contamination in the maps. 
Since we seek to derive the most stringent upper limit 
on r that Spider could set given the assumptions in this 
paper, the input cosmological model assumes r = 0.001 
(and r = 0.09). Farhang et al. (submitted) address the 
question of how well "large" -r £?-modes can be detected. 



With foregrounds 
No foregrounds 




0.04 



Fig. 8. — Marginalized one-dimensional r-likelihood curves for 
Spider and Planck HFI. In the "with foreground" case, foreground- 
subtraction residuals are included as contaminants. The "no fore- 
ground" case assumes no foreground contamination. In the absence 
of foregrounds, the first flight of Spider (Spider 1) achieves the 
3cr limit r < 0.03, while the two flights combined (Spider 2) lead 
to r < 0.02 (99% CL). In the "with foreground" case, two nights 
are needed to reach the 3cr limit r < 0.03. The details of the two 
flights are given in Table [4] 

We find that, with two flights, Spider has the abil- 
ity to constrain the value of the tensor-to-scalar ratio 
r to below 0.03 (99% CL). This number improves to 
0.02 (99% CL) in the limit of no foreground contamina- 
tion. It is worth noting that the constraint quoted above 
for the "with foreground" case is based on the foreground 
removal technique described earlier in this section, which 
assumes that the frequency dependence of the polarized 
dust emission is well described by a single power law. 
Although this is in line with t he current state-of-th e-art 
for _B-mode forecasts (see, e.g., Dunkley et al.|2009| , the 
detailed question of how deviations trom this description 
will be handled has yet to be studied. Any such deviation 
will likely impact the quoted constraints. Relatedly, it is 
possible that a different frequency dependence of the po- 
larized dust emission might require a different frequency 
coverage in order to reach the best possible constraints 
on r. Should this become apparent after the first Spi- 
der flight, or through improvements in our knowledge 
of foregrounds before Spider flies, the flexibility of Spi- 
der's focal plane composition discussed in Section 2.1 
ensures that we will be able to adopt the best frequency 



coverage to reach our science goals. 

6. SUMMARY & CONCLUSIONS 
Our main results are as follows. 

1. The Spider instrument was designed with an em- 
phasis on reaching low levels of polarized system- 
atic error at multipoles close to the I ~ 80 ac ousti c 
peak in the B-mode power spectrum (Section |2.1| ). 
Taking into account all known sources of system- 
atic error, the residual power spectrum is less than 
37% of the B-mode signal for r = 0.03 at I = 100 
(Section 2.2 1. As we refine our systematics budget, 
we expect this number to go down for three rea- 
sons: (1) more stringent design targets than those 
shown in Table [2] have already been reached; (2) re- 
cent updates to our flight plan and to our scanning 
strategy are expected to reduce the overall level of 
systematic error by providing better observing con- 
ditions; and (3) all residuals in Table [2] were added 
in quadrature, which leads to a conservative esti- 
mate of the global systematic error. Simulations 
taking these three points into account are under- 
way. The reduction in systematic error that we ex- 
pect these improved simulations to show will likely 
be partially canceled by the inclusion of systematic 
effects that are yet to be characterized, in partic- 
ular in so far as they relate to thermal stability of 
the focal plane and of the optical sleeve. 

2. Polarized sidelobes should not hinder Spider's 
ability to detect or constrain £?-modes. Using the 
polarized sidelobe profile measured in Bicep, we 
showed that the corresponding contamination at 
£ = 100 is over an order of magnitude fain ter than 
the r = 0.03 £?-mode signal (Section |2.3[ ). As il- 
lustrated in Figure [2] improved baffling will lower 
this number by over an order of magnitude. At the 
largest scales [l ~ 10), it will lead to a polarized 
sidelobe signal of amplitude comparable to that of 
the r = 0.03 £?-modes. 



3. Spider's observing strategy (Section 4.1) was de- 
signed to minimize the T/Q/U covariance, thereby 
minimizing the uncertainty in the reconstruction of 
the Stokes parameters. Figure [5] shows that, taking 
all detectors into account, the excess variance in- 
troduced by the scanning strategy is less than 1% 
over 96% of the sky observed during a 4-day period. 

4. The lack of atmospheric contamination at strato- 
spheric altitude and large sky coverage enabled by 
Spider's small apertures provide a relatively broad 
dynamic range in multipole, as shown b y th e filter 
transfer function (Ft) study in Section 4.2 This 
opens the possibility of moving beyond detection 
to a spectral characterization of the B-mode power 
spectrum, an important consistency check on the 
cosmological origin of any detected signal. Further- 
more, the fidelity of the measurement results in a 
significantly higher final sensitivity than would be 
achieved by an experiment of the same NET using 
a less efficient observation platform. 
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5. The frequency distribution as a function of Spider 
flight summarized in Table [4] ensures high signal-to- 
noise measurements of the microwave sky at low- 
foreground frequencies, while the polarized emis- 
sion from interstellar dust is mapped with high 
sensitivity in a band where it overwhelmingly dom- 
inates. With this distribution, we showed in Sec- 
tion [5] that Spider will provide the stringent con- 
straint r < 0.03 (99% CL) after both flights, even 
when accounting for residual foreground contami- 
nation. This number is derived under the assump- 
tion that the Planck 217 GHz data will be avail- 
able during the Spider analysis. When foreground 
contamination is ignored, Spider's constraint im- 
proves to r < 0.02 (99% CL). 

6. Among Galactic foregrounds, we expect the polar- 
ized emission from interstellar dust to be the only 
signal as bright or brighter than the r — .03 B- 
modes at Spider frequencies (Section 3.1). We 
model this emiss ion following the procedure de- 
tailed in Section |3.2| and show that the Spider 
observing region includes the cleanest 2% of the sky 
accessible from a McMurdo LDB flight. The aver- 
age amplitude of the polarized Galactic dust emis- 
sion in the Spider field is comparable to that in a 
typical / s ky = 2% patch in the "Southern Hole;" 
many / s k y = 2% patches in the Spider region 
exhibit significantly less polarized dust emission. 
This study suggests that without component sep- 
aration, degree-scale polarized dust emission will 
limit the constraints of any experiment at or above 



the level of r ~ 0.03, even in the portions of the 
southern sky most free of Galactic dust emission. 

Although more work is needed to fully characterize the 
Spider instrument and its nontrivial interaction with 
the microwave sky, a simple Galactic foreground situ- 
ation could allow Spider to characterize a cosmologi- 
cal -B-mode signal whose intensity is close to the lowest 
compatible with the simplest viable class of inflationary 
models. General studies of the impact of a more complex 
foreground situation have yet to be performed, and are 
of obvious interest in the context of Spider. 
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